Emergent thermodynamics in a quenched quantum many-body system 
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We study the statistics of the work done, the fluctuation relations and the irreversible entropy production in a 
quantum many-body system subject to the sudden quench of a control parameter. By treating the quench as a 
thermodynamic transformation we show that the emergence of irreversibility in the nonequilibrium dynamics of 
closed many-body quantum systems can be accurately characterized. We demonstrate our ideas by considering a 
transverse quantum Ising model that is taken out of equilibrium by the instantaneous switching of the transverse 
field. 



Introduction- In the past decade or so, there has been a revival 
of interest in the study of nonequilibrium dynamics in closed 
quantum systems. Mainly, this is due to a series of spectacular 
experiments using ultracold atoms, where the high degree of 
isolation and long coherence times permit the study of dynam- 
ics over long timescales (TJ. These experiments have raised 
a number of important theoretical issues including the rela- 
tionship between thermalisation and integrability [2| and the 
universality of defect generation in the adiabatic crossing of a 
critical point |3 1. A common way to take a many-body system 
out of equilibrium is by an abrupt change of a local or global 
parameter of the Hamiltonian, this is commonly referred to 
as a 'sudden quench'. Following a quench the dynamical re- 
sponse of the system can be probed by studying, for example, 
the dynamical correlation functions [4], change in the diago- 
nal entropy [5 1 or the statistics of work done |6|. 

Over a similar period of time, there has also been a great 
deal of interest in the statistical mechanics community sur- 
rounding the discovery of the nonequilibrium fluctuation rela- 
tions (see e.g. Ref. [7| for a review). Essentially, the fluctua- 
tions relations encode the full non-linear response of a system 
to a time dependent change of a Hamiltonian parameter. In 
particular, they make a definitive statement regarding the irre- 
versible entropy production of a system following a thermody- 
namic transformation and, as such, allow us to understand the 
emergence of thermodynamic behaviour in systems where the 
microscopic laws are inherently reversible. Given the current 
experimental interest in the nonequilibrium dynamics of ultra- 
cold atomic systems and the recent developments in statistical 
mechanics, it is natural to study the quench dynamics of quan- 
tum many-body systems in this new thermodynamical formu- 
lation. In this work we use the transverse Ising model 1 8 1 to 
provide an exact analysis of the Tasaki-Crooks and Jarzynski 
fluctuation relations in a quenched many-body system. Fur- 
thermore, we compute the average irreversible production and 
show that emergence of thermodynamics provides an elegant 
interpretation of the essential physics. 



Nonequilibrium quantum thermodynamics- We begin by re- 
viewing some key concepts of microscopic thermodynamics, 
allowing us to define the formalism that is used in the rest of 
our study. 

One of the fundamental goals of quantum thermodynam- 
ics is to understand how thermodynamical laws emerge from 
the underlying quantum mechanics of individual particles f9|- 
In this spirit, we consider a dynamical system described by 
a Hamiltonian H(X(t)) that depends on an external work pa- 
rameter A(i), i.e. an externally controlled parameter whose 
value determines the equilibrium configuration of the system. 
The system is prepared by allowing it to equilibrate with a 
heat reservoir at inverse temperature /3 for a fixed value of the 
work parameter A(t < 0) = Ao- The initial state of the system 
is thus the Gibbs state Pg(Ao), where 

e -pH(X) 

Pg{x) : = ~z(xj-> 

and the partition function Z(X) = Tr [e~ /3H(A) ]. At t = 
the resevoir-system coupling is removed and a protocol is per- 
formed on the system taking the work parameter from its ini- 
tial value Ao to a final value A T at a later time t = r. The initial 
and final Hamiltonians connected by the protocol Ao — > A T 
have the spectral decompositions ff(Ao) = J2n e n(Xo)\n) (n\ 
and H(X T ) = 5Z m eJn(A T )|m)(m|, respectively, where \n) 
(|m)) is the n th (m th ) eigenstate of the initial (final) Hamilto- 
nian with eigenvalue e„ (e' m )- 

The definition of the work done on the system If as a 
consequence of the protocol requires two projective measure- 
ments; The first projects onto the eigenbasis of the initial 
Hamiltonian H (Ao) at t = 0, with the system in thermal equi- 
librium. The system then evolves under the unitary dynamics 
U(t; 0) generated by the protocol Ao — > X T before the second 
measurement projects onto the eigenbasis of the final Hamil- 
tonian H(X T ). The probability of obtaining e n for the first 
measurement outcome followed by e' m for the second mea- 
surement is then P^p T m \ n = e~^ e ™ \{ti\U(t, 0)|m) | 2 /Z(A ). 
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Accordingly, the work distribution is defined as ITOl 
P ¥ (W) := Y^Pn P T m\n 6 l W ~ (4 - *»)] • 



(1) 



Eq. ([T| therefore encodes the fluctuations in the work that arise 
from thermal statistics (p°) and from quantum measurement 
statistics (p^ n ) over many identical realisations of the proto- 
col. For our purposes, it is convenient to define the character- 
istic function of the work distribution as the Fourier transform 
QfEq.QEl 



Xf(u,t)= I dWe mW Pp(W), 



Tr 



U\t, 0)e tuH{K) U(T,0)e~ luH(Xo) p c (Xo) 



(2) 



The convenience of x F (u,r) is evident when consid- 
ering the well-known Tasaki-Crooks fluctuation relation 
P F (W)/P B (-W) = e /3(^-AF) fiol [121. This states that 
the ratio between the forward work distribution Pp(W), intro- 
duced above, and the backward work distribution Pb(— W), 
obtained from the protocol A T — > A in which the system is 
initialized at t = in the Gibbs state pg(A t ) and evolves ac- 
cording to W (r, 0), is related to the difference in the equilib- 
rium free-energy of the system AF. Following Ref. Ifl3l the 
Tasaki-Crooks relation is written in terms of the characteristic 
function as 



XF(u,T) = Z(\ T ) 

Xb(v,t) Z(A ) : 



(3) 



where we have introduced the backward characteristic func- 
tion xb{v) = J dW e lvW Pb{W). Moreover, the Jarzynski 
equality lfT4l is easily obtained from Eq. |2]) by introducing 
the parameter u — i/3, giving 



Xf(*/3,t) 



-pW\ 



Z(X r ) 

Z(X ) 



-0AF 



(4) 



where in obtaining the last equality we have used the relation 
AF = -(l//3)ln(Z(A r )/Z(A )). Both the Tasaki-Crooks 
and Jarzynski fluctuation relations are statements regarding 
the symmetry of fluctuations in work during thermodynamic 
transformations of microscopic systems. Remarkably, these 
symmetries are solely determined by the equilibrium state 
quantity AF regardless of how far the system is driven from 
equilibrium. For a recent information-theoretic interpretation 
of the fluctuation relations see Ref. ifTSl . 
Irreversible entropy production- For finite systems, the statis- 
tical nature of work Eq. ([TJ requires the second law of thermo- 
dynamics to be revised to the form (W) > AF, with equality 
being reached for a quasistatic process. For all non-ideal pro- 
cesses, the deficit between average work (W) and the varia- 
tion in free energy can be accounted for by the ad hoc intro- 
duction of the average irreversible work, 

(W) = (Wfa) + AF. 



For a closed quantum system, the heat transfer into the sys- 
tem (Q) = and the sole contribution to the change in en- 
tropy is the average irreversible entropy production A5; 1T . In 
Ref. |fT6ll it is shown that for an initial Gibbs state pg(Xq) 
undergoing unitary evolution generated by a time dependent 
Hamiltonian H(X(t)), the average irreversible entropy pro- 
duction is given by the relative entropy between the instanta- 
neous state of the system p(t) = U(t, Q)pc(Xo)U^ (t, 0) and a 
hypothetical Gibbs state at that time, 



(AS llT ) = S (p(t)\\ Pc (X t )) 

= Tr \p(t) logp(i)] - Tr [p(t) log p G (A(i))] 



(5) 



In the case of a sudden quench, in which the work parame- 
ter A(t) is suddenly switched between some initial and final 
value, Eq. |5) becomes 

(AS m ) = S(p c (X )\\p G (X T )) = 0({W) - AF) - p{W m ). 

(6) 

This expression for the average irreversible entropy produc- 
tion induced by a sudden quench in a closed quantum system 
was first discovered by Donald in Ref. ifTTl within a different 
context. 

Transverse Ising model- We now apply the framework of 
nonequilibrium statistical mechanics outlined above to the 
nonequilibrium transformation of a thermal quantum spin 
chain. In particular, we analyse the sudden quench of the 
transverse field in the quantum Ising model. For a discussion 
of this model in the zero temperature limit see Refs. l6l [T8l . 

We consider a one-dimensional ring of TV spin- 1/2 particles 
that interact with their nearest-neighbours via ferromagnetic 
coupling along the z-axis and with an external field applied 
along the x-axis. The Hamiltonian is 



N 



(7) 



where A is a dimensionless parameter measuring the strength 
of the external field with respect to the spin-spin coupling, cr™ 
(a = x, y, z) is the spin-1/2 Pauli operator acting at the j th 
spin and periodic boundary conditions are imposed by requir- 
ing that cr^r +1 = erf . The Hamiltonian Eq. |7| is diagonalised 
by decomposing the Hilbert space into orthogonal parity sub- 
spaces and following the procedure outlined in the supporting 
material. In this way, considering the positive parity subspace 
only, the initial Hamiltonian with A = Ao is written E2l 



H(X ) 



lllk 



(8) 



where jk,lt are fermionic creation and annihilation op- 
erators labelled by the members of the set K + = 
{±7r(2n — 1)/N : n = 1, ...N/2} of positive parity subspace 
pseudomomenta. Proceeding as earlier, the system is prepared 
in the Gibbs state pg(A ) = e~ /3ff(Ao) /2(Ao) with inverse 
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temperature (3 and associated partition function 
Z(X )=2 N J] cosh 2 ^ £fe(Ao) 



k£K+ 
k>0 



The protocol constitutes the instantaneous switching of the 
transverse field to the final value A = A T , giving the final 
Hamiltonian 



H(X T ) = £ < X t) 



k£K+ 



(9) 



Note that the differing values of the transverse field in Eqs. <(8j 
and |9]l require diagonalising transformations that are quan- 
titatively different. Consequently the post-quench fermionic 
operators {7fe} differ from their pre-quench counterparts 
{jk}, though the allowed values of the psuedomomenta are 
identical in both cases. In the case of a sudden quench the 
characteristic function Eq. |2) takes the simplified form 



Xf(u) 



Tr 



JuH(X T ) -iuH{\ ) + 



(10) 



Using Eqs. ([8) and (|9j, the trace in Eq. ( [10] ) is taken over the 
eigenstates of the initial Hamiltonian {|n&, n^ k )} to give 



Xf(w) 



Z(A 



n e ■ 

fe>0 



-(ra+/3)efc(Ao)(rafc+ri_ fc -l) 



(n fe ,n_ i; |e i " e '' (v) (^+^-^-*- 1 )|n i: ,n_ fe ). 



The matrix elements can be evaluated explicitly to give an 
analytic form of the forward characteristic function. Hence, 

fe>0 

+ e -(i«+«e fe (A ) (c+( Ui A T ) + ^(u,A T )) +2. 

(11) 



Here 



we 



have 



introduced 



the 



quantities 



Cf(u,A) = cos 2 (A fc /2)e ± " 1£fc ( A ) and Sjf(u.A) = 
sin 2 (A /£ /2)e ±iu£fc ( A \ where A fc = <t> k - fc is the differ- 
ence in the pre- and post-quench Bogolyubov angles (see 
supporting information). 

Verification of the fluctuation theorems- The verification of 
the Tasaki-Crooks relation Eq. <(3j requires an expression for 
the backward characteristic function. This is easily obtained 
using a procedure similar to that described above for the 
forward characteristic function Eq. ( fTT| under the mapping 
A O A r e fc (A ) o e fe (A T ), A fc -> -A fc . The Tasaki- 
Crooks relation follows from Xb( v ) by introducing the com- 
plex parameter v = —u + Noting that (—u + i/3, A) = 
C%(<u,\)e*P"-W and S± {-u + i/3, A) = Sg='(«, A)e=F^W, 
it is straightforward to show that 



Z(X T ) 



k£K + 
k>0 



+ e -(i«+«e*(Ao) (c+( u ,A T ) + 5fe (u,A T )) +2. 



The ratio of the forward and backward characteristic functions 
is thus 

Xf(u) _ Z{X T ) 
Xb{v) Z(X )' 

which is equivalent to the Crooks relation Eq. (BJ. Further, the 
Jarzynski equality Eq. Q follows from the forward character- 
istic function Eq. ( fTT| ) by introducing the complex argument 

u = i/3, 

X*W)=zfx) 1] (2 + 2cosh(/3e fc (A r )) 

° k£K+ 
k>0 



Z(X ) 



n c ° sh2 



fee if + 
k>0 



f3tk{Xr) 



Z(X T ) 

Z(X Q )' 



To our knowledge this is the first analytic demonstration of 
the fluctuation relations in a non-trivial quantum many-body 
system incorporating a critical point. 

Emergent thermodynamics- The general form of the forward 
characteristic function following a sudden quench Eq.(fT0|i ad- 
mits a simple expression for the average work, 

(W) = Tr [H(X T )p G (X )] - Tr [H(X )p G (X )] . (12) 

Using the approach presented in the supporting material, the 
evaluation of Eq. ( fT2] i leads to the following closed analytic 
form for the average work done on the spin system 



(W)= MA )-e fe (A T )cos(A fe ))tanh 



k£K+ 
k>0 



(3e k (X ) 



= 2(A -A r ) 



V] cos (<p k ) tanh | ^ £fc ( A °) j _ 
\ z / 



keK 
k>0 



This in turn allows the calculation of the irreversible entropy 
production Eq. |6} for arbitrary temperature, spin number and 
quench amplitude, 



(A Sin-) = P(W) 



E 

keK+ 
k>0 



hi 



cosh 2 C8e fc (A T )/2) 
cosh 2 (/?e fe (A )/2) ' 



This quantity and its thermodynamic interpretation encap- 
sulate the physics of the quench problem in a closed critical 
system in a remarkably simple way. Fig. [T] shows the irre- 
versible entropy production due to a series of sudden quenches 
with amplitude |A r — A | = 0.01. The left figure shows the 
quantity for spin chains of several sizes and low temperature 
f3 = 100. The interpretation of the behaviour is straightfor- 
ward; As the size of the system increases the energy gap be- 
tween the ground and first excited state at the critical point 
begins to close. Work is performed to drive the system across 
the critical region and, due to the vanishing energy gap, it 
becomes increasingly difficult to do so without exciting the 
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A. ; =X +0.01; p=100 X =^+0-01- N=10,000 




\ ' \ 

FIG. 1: Left: The irreversible entropy production for a series of 
quenches with amplitude |A T — Ao = 0.01 at /3 = 100 for sev- 
eral different ring sizes. The initial value of the transverse field Ao 
is shown on the a>axis. Right: The irreversible entropy production 
for the same series of quenches in a ring of N = 10, 000 spins at 
various temperatures. 



system, thereby dissipating work. This leads to the produc- 
tion of irreversible entropy and the emergence of intrinsic ir- 
reversibility in the critical region. 

The figure on the right shows the irreversible en- 
tropy production in a chain of N = 10, 000 spins at 
various temperatures. As expected, the signature of 
quantum criticality decreases at higher temperatures with 
the emergence of thermal fluctuations. The source of 
irreversibility is elucidated by manipulating the Tasaki- 
Crooks fluctuation relation to obtain the expression 
(AStr) = J AW P F (W)log(P F (W)/P B (-W)) = 
K{P V {W)\\P B {W)), where IC {P F {W)\\P B (W)) is the 
Kullback-Liebler relative entropy, measuring the distance 
between two probability distributions. Intuitively, this 
expression attributes the uncertainty in distinguishing the ex- 
perimental data contained in the forward and backward work 
distributions to the degree of irreversible entropy production. 
According to this powerful interpretation, quantum criticality 
has the effect of setting the thermodynamic arrow of time 
as the irreversible entropy production grows with decreasing 
temperature. 

As a final remark, we note that the forward characteristic 
function Eq. OJ) has a similar form to the Loschmidt echo 
lfl9l . This has been shown in previous work to be a good in- 
dicator of phase transitions E51 and could be experimentally 
measured using Raman interferometry EH . 
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Appendix A: Diagonalization of the transverse Ising model 

The quantum Ising model in a transverse field describes 
a lattice of spin- 1/2 particles that interact with their nearest- 
neighbours via ferromagnetic coupling along the z-axis and 
with an external field applied along the x-axis. For a spatially 
homogeneous one-dimensional lattice of N spins in a uniform 
field, the Hamiltonian is 

N 
3=1 

where A is a dimensionless parameter that measures the 
strength of the external field and the Pauli spin- 1/2 operators 
are defined with periodic boundary conditions cr^ r+1 = a" 
(a = x, y, z). Under the canonical transformation <rj — > 
erj, crj — > — crj Vj the Hamiltonian Eq. ( |S-l| l becomes 

N 

H = -J2^j+°M+i, ( S "2) 

3=1 

The spin operators are mapped to a spinless fermionic opera- 
tors by a Jordan-Wigner transformation, thus 

c j = 2 n^K+^P' 

i=i 

i=i 

Here, the operators Cj (cj) annihilate (create) a Jordan- Wigner 
fermion at the j th lattice site and obey the usual fermionic anti- 
commutation relations. This in turn allows the definition of 
the parity operator 

N 

n : = n i 1 2c h) . 

j=i 

which measures whether the number of fermions in the chain 
is even (n = 1) or odd (n = —1). Following the Jordan- 
Wigner transformation the Hamiltonian Eq. ( |S-2| > factorises 
into two orthogonal parity subspaces, 

H = P+H+P+ + P~H~P~. 

Here, P = (1/2) (1 ± n) are the projectors onto the even 
(+) and odd (-) parity subspaces and 

N 

H ± = -X-Y / (2Ac] Cj - (c]c j+1 + c J+lCj + h.c.)) , 

(S-3) 

are the even and odd parity subspace contributions to the 
Hamiltonian. The even and odd parity subspace Hamiltoni- 
ans are identical with the exception that in H + we impose 
the boundary condition cjv + i = — ci and in H~ we impose 
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cn+i = c\. Note that the parity of the chain is conserved 
, II] = 0); the Hamiltonian in Eq. ( |S-3| l does not mix the 
parity subspaces. Initialising the system in a state with zero 
projection onto, say, the odd subspace then restricts the dy- 
namics to the even subspace only. For an initial Gibbs state, 
the system is a mixture of positive and negative parity states 
and both subspaces must be accounted for. Despite this we 
restrict our attention to the even parity subspace only. This 
treatment becomes exact in the thermodynamic limit where 
boundary effects become negligible. This is the correct limit 
in which to discuss phase transitions, however a discussion of 
the fluctuation theorems is more suited to a finite chain. The 
analysis of a chain of arbitrary length in which both parity sub- 
spaces are accounted follows from a straightforward extension 
of what is presented here at the expense of more cumbersome 
expressions and provides little extra insight. Note that, in the 
main text, the fact that we consider the even subspace only is 
denoted by the summation over the set of positive parity pseu- 
domomenta k G K + in all relevant expressions. In the main 
text the '+' superscript that is used here to explicitly distin- 
guish between the positive and negative parity subspaces is 
dropped for the sake of brevity. 

The diagonalisation of the Hamiltonian is completed by 
application of a Fourier transformation followed by a Bo- 
golyubov transformation, which factorizes over the spaces 
with different pseudomomentum k. For the positive parity 
contribution to the Hamiltonian the Fourier transformation is 
defined as 



t/4 



N 



E 

k£K+ 



c k e 



where the values of pseudomomentum are 

K+ = {k = ±^{2n-l), n = l,...,y 



After this transformation the Hamiltonian Eq. ( |S-3| l takes the 
form 

H+ = 2(A-cos(fe))4c fe +sin(fc)(c|.cL fe +c_ fc c fe )-A. 

Note that all the terms preserve pseudomomentum so that the 
remaining step of the diagonalization can be performed within 
each subspace with assigned value of ±k. The last step is the 
Bogolyubov transformation 



c±k = l±k cos I — I T %k sm ( Y 



(S-4) 



where 



cos(</> fc ) 
sin( < /)/ c ) 



A — cos(fc) 



y / sin 2 (fc) + [A-cos(fc)] 2 ' 
sin(fc) 

^/sin 2 (fc) + [A-cos(A:)] 2 



(S-5) 



With this, the Hamiltonian can be written in the form 

h + = E ^ (iH - J' 



keK+ 

with the dispersion relation 



e k = 2yW(» + [A - cos(fc)] 2 . 

Note that e k = e-k > and that the total spectrum is sym- 
metric with respect to the zero of energy. 



Appendix B: Connecting the initial and final Hamiltonians 

To evaluate the characteristic function explicitly, the eigen- 
states of the initial Hamiltonian H + (Xq) must be written in 
terms of the eigenstates of the final Hamiltonian H + (\ T ). In- 
verting Eq. \S-A\ and its hermitian conjugate it is possible to 
relate the sets of pre- and post-quench Bogolyubov operators. 
Hence, 



,'A fe \ t • f A * 



1-k = 1-k cos 



A, 



Here Afc = <f>k — 4>k and the expressions for the pre- and post- 
quench Bogolyubov angles, 4>k an d </>fc> nave tne form given 
in Eq.( |S-5[ ) with A = A and A T respectively. Using this, the 
vacuum states in the two representations are related by 



|0 fe ,0_ 



Afc 
2 



Afc 
2 



7fc7lfcj |6fc,6_fc). 

(S-6) 

The expressions for higher energy eigenstates \n k ,n_k) are 
then obtained by applying the appropriate creation operators 
to Eq. dS^6). 



Appendix C: Calculation of the average work 

The calculation for the average work done on a quenched 
transverse Ising model takes advantage of the factorization of 
_ff + (A ) and H + (X ) into blocks of paired pseudomomenta 
with labels ±k. We start by writing the density matrix of the 
system as p = ®fc>o P±k> so that 



(W) = Tr 



£(i/+fc(A T )-i^ fc (Ao))(g)p ± fc 



(S-7) 



where H^_ k (\ T ) = efc(A r )(7^7 fe + rLfcT-fc - 1) and 
i^fc(A ) = £fc(A )(7fe t 7fc + 7-fcVfc - 1). With a little 
effort, Eq. (|S-7[> can be rewritten as 



k>0 



(H+ k (X T )-H+ k (X ))l[a ±k 



k'>Q 



, (S-8) 
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where, 



o±fc = ^ \n k ,n- k )(n k ,n_ k \ 



n± k =0,l 



; -^e fe (Ao)(nfe+n_fe-l) 

4cosh 2 (/3e fe (A )/2) ' 



Noting that Tr [nfc'>o °"±fc'] = 1* Eq.( |S-8| l reduces to the 
form (WO =E fc >o( M/ fc)' with 

(W*) = Tr[(H+ k (X T ) - H+ k (X ))a ±k ]. 

In order to calculate the trace we note that we need only 
keep the terms of H± k (\ T ) that are diagonal in the basis of 

[i^ fc (A T )] diag = e fe (A T ) cos(A fe )( 7fc t 7fc + 7 _ fc t 7 _ fc - l). 
With this, (Wfe) takes form 

(W k ) = (cos(A fe )e fc (A r ) - £fc(A )) Tr [(n k + n_ fe - l)a± k ] , 

which leads straightforwardly to the expression for the aver- 
age work in the main text. 
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